THE BULLETIN OF 


Mathematical 
BIOPHYSICS 


HE UNIVERSITY OF CHICAGO PRESS - CHICAGO - ILLINOIS 


-OLUME 3 
}UMBER 3 
EPT. 1941 


THE BULLETIN OF 
Mathematical Biophysics 


ErDel TED BeYs @N 2 R BA GSEH E V S-KsY 


The Bulletin is devoted to publications of research 
in Mathematical Biophysics, as described on the 
inside back cover. 


THE BULLETIN is published by the University of Chicago at 
the University of Chicago Press, 5750 Ellis Avenue, Chicago, IlIli- 
nois, quarterly, in March, June, September, December. {[The sub- 
scription price is $2.50 per year, the price of single copies is 75 
cents. Orders for service of less than a full year will be charged at 
the single-copy rate. [Patrons are requested to make all remittances 
payable to The University of Chicago Press in postal or express 
money orders or bank drafts. 


THE FOLLOWING are authorized agents: 


For the British Empire, except North America, India, and 
Australasia: The Cambridge University Press, Bentley House, 
200 Euston Road, London, N.W. 1. Prices of yearly subscriptions 
and of single copies may be had on application. 

For Japan: The Maruzen Company, Ltd., Tokyo. 

For China: The Commercial Press, Ltd., 211 Honan Road, 
Shanghai. 


CLAIMS FOR MISSING NUMBERS should be made within the month 
following the regular month of publication. The publishers expect 
to supply missing numbers free only when losses have been sus- 
tained in transit, and when the reserve stock will permit. 

BUSINESS CORRESPONDENCE should be addressed to The Uni- 
versity of Chicago Press, Chicago, III. 

COMMUNICATIONS FOR THE EDITOR and manuscripts should be 
addressed to N. RASHEVSKY, Editorial Office of the Bulletin of 
Mathematical Biophysics, 5822 Drexel Avenue, Chicago, III. 


—_—_— eee 
PRINTED BY THE DENTAN PRINTING CO., COLORADO SPRINGS, COLORADO 


BULLETIN OF 


Mathematical Biophysics 


CONTENTS 


ELECTRICAL CHARGES AND POTENTIALS IN CELLS 
RESULTING FROM METABOLISM OF ELEC- 
PROD YER = “soe yy Soe ey se a es TD 
ROBERT R. WILLIAMSON 


NOTE ON THE MATHEMATICAL BIOPHYSICS OF TEM- 
PORAL SEQUENCES OF STIMULI - - - - 89 
N. RASHEVSKY 


A NOTE ON THE NATURE OF CORRELATIONS BETWEEN 
DIFFERENT CHARACTERISTICS OF ORGAN- 
DMs eS cm ea A er we 98 

N. RASHEVSKY 


MATHEMATICAL BIOPHYSICS OF THE GALVANIC SKIN 
RECON TMM ees gee ey 97 
CLYDE H. COOMBS 


A THEORY OF STEADY-STATE ACTIVITY IN NERVE- 
FIBER NETWORKS: II. THE SIMPLE CIRCUIT - 105 
ALSTON S. HOUSEHOLDER 


VOLUME THREE SEPTEMBER 1941 NUMBER THREE 


is i 


= 4 = 


2 TRAE 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 3, 1941 


ELECTRICAL CHARGES AND POTENTIALS IN CELLS 
RESULTING FROM METABOLISM OF ELECTROLYTES 


ROBERT R. WILLIAMSON 
THE UNIVERSITY OF CHICAGO 


An approximate solution for the relation of the charges and po- 
tentials in a spherical cell to the production of electrolytes in the cell 
and their diffusion resistances is derived. The potential obtained by 
introducing reasonable values of the constants is of the proper order 
of magnitude. The equations are applied to a respiratory chain and 
the relation between oxygen consumption, glycolytic coefficient, and 
potentials is determined. Available experimental data is compared 
with the theory. 


In the following derivations we will use these notations. 
D.= diffusion coefficient of the positive ion in cm? sec, 
D_.= diffusion coefficient of the negative ion in cm? sec", 


h.—= membrane permeability of the positive ion in cm 
sec, 


.. = membrane permeability of the negative ion in cm 
sec —, 


Q = rate of production of ions in mol cm- sec, 

C = internal concentration in mol cm, 

C, = external concentration in mol cm~, 

1) = radius inside membrane in cm, 

6 = membrane thickness in cm, 

T = absolute temperature in deg A., 

f = electrical force on an ion due to field of other ions. 
k = Boltzmann constant = 1.4-10-** erg deg, 

e= charge on an electron = 4.8-10-” abs. E. S. units, 
p= charge density, | 
K = dielectric constant, 
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4 wt Pee N e 
Kkt 
ly. (Williamson 1939), 


N = Avogadro’s number = 6.06-10?? molecules mol", 


A= = 4.7.10", with values chosen previous- 


J = transport in mol cm” sec’, 

B= mobility of an ion in cm? erg” t", 

A, = diffusion resistance of the positive ion in sec, 
See Lo 


= De take 


A_= diffusion resistance of the negative ion in sec, 


Pe: "potential difference between 7, and 72 in abs. E. 8. 


volts. 


Denote by a bar average, by a + or — subscript the ions to which 
the quantity refers, and by a subscript 1, values at the membrane, by 
subscript 0, external values, except 7, , radius of cell. 

Assume a cell of two phases, spherically symmetric, one extend- 
ing from 0 to 7, and the other constituting a shell of thickness 6 corre- 
sponding to a membrane. Assume a non-electrolyte to be metabolized 
producing an electrolyte. If the diffusion resistances are different for 
the ions of opposite charge, there will be a difference in concentration 
which will result in a net charge. If we assume the distribution of 
the ions to be linear, the charge density, given at a point by 


o=NeA, where A=C,—-C., (1) 
will also be a linear function of r. 
Let 
p—ar+b,where0 <7 <7, (inner phase), 
and oo. (2) 
p= =a,r + b, wherer, <7 < 1% i 6 (outer phase). 


The force, resulting from the charge, on a unit charge ata a Eon ris 
given by: (Joos, 1934) °' 


F=7-, , where E = An | p(n) rar. % ye 


Fae 
% ) 
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The potential difference between any two points 7, , 7. is 


PD: 


i f "Fdr. (4) 


In the inner phase 


V2 


Ser Gs (oe : 

P.D. = r{=i, r? p(r) dr, 
(5) 

4nf a A < bo ae 
me ae 
and if 7, =0,andr,=%, 
To me ig 

ay ae a gr (6) 


3K 


0 
if we let A= A at r,/2. 


In the membrane, equation (5) applies, but p(r) = ar + b or 
ar + b, according to equations (2). Therefore, 


r a-— @ b= <b agr*-b4r* 
72. 9(r)dr = Trt + 73 + = Z ; 
J o(7) TD sore tts 
(7) 
To+d 10 
PD, = P.D. =x Bar,” ae Abr. | ’ 


when second and higher powers of 6/7, are neglected. 

A is a function of r. By the standard reasoning used in the ap- 
proximation method (Rashevsky, 1940), we consider that C has ap- 
proximated the values C when r = r,/2. Similarly we define C, as 
the value of C for r—vr,. Then let 


AC. Ce, 
Ape Ce 3G fe 
and let Du. = De- = ~. 
Then from equations (1), (2), and (8) 
$= mal (A — A,), and 
E (9) 
b=Ne(2A—A,). 
Introducing these values into equation (7) we have 
Pp nee aad. (10) 


m 
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Now we must evaluate 4, A,, from the diffusion equations. In 


the stationary state the material balance equations are 


Keds =I: se (11) 


The equations for the transport may all be expressed in the form 
(Reiner, 1937, Williamson, 1939) 


J— 1) e740 Gee Bice 


where (Gyemant, 1925, p. 55) B= D/kT. 

If we take an average gradient and force on the ions (N. Rashev- 
sky, 1940), and introduce the subscript notation we obtain from the 
above equation the following four steady state equations for trans- 
port of respectively, positive ion in the inner phase, positive ion in 
the outer phase, negative ion in the inner phase, and negative ion in 
the outer phase, after slight rearrangements. 


0. Cn eee (C.—C_)C,, (12) 

oo CeCe eee (C.— C (C2 C,.) ae 
and 

= C.-C. — es (CocaC.) Ca, (14) 

we yeti ies oe (COSC CST 


Adding equations (12) and (13), and (14) and (15), we have 
ra FE Oe x 
.Q=6,— Cn + (G,-E) [é. +5 (C+ Cu) |, 
bs A ik ode = 
AQ=C_—C,.= = (C,— C_) |C. + RS (Ci. + Cc.) | : 
Or since A © 10%, 6/7, © 10° or less 
Ass. ae 
MO) (C= CIC, LO== 
Cn Qa (C.2— C2). 


If we introduce the notation 
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C=C AS C=C +A. , expand, and drop powers of A, or A- 
higher than the first, we have 


—___ Q=4,- 4.=A. (16) 


Collecting C,,, Ci. in equations (13) and (15) we have 


0) = Be 
(ee Wa yet © [i -% 4], 
27> 


3h 2ro 
AéA Qro AéA 
(1- = a) Oe = ay + Oo [t+ i |. 


If we divide by the coefficient of C,, apply the binomial expansion 
and consider only the first two terms of the expansion, we have 


ee Se [tS | ofA], 


3h. 27 27 
Q7ro AéA AéA ia 
C= 5" [1+ fs | + [1+ cae 


or 


A, =— 


Qro fly —h., ASA bth.) _2C,AdA4 
3 UE eae a ‘oe: 


Choosing the following plausible values, (Rashevsky, 1940; R. Hober, 
1936) 


h,=10- em? sec?;  h_—=10"cm?sec?; 6=5-107 cm; 
Ao secs, AS 200 gees” A— 100 sec; K = 80; 


Co=10- mol em;. Q—= 107 mol cm= sec“; 


we see that ; , 
ihe — Ri AO eR he 2A6A 
s = Se Oey 
hl or, Wah. ares 
and we may write this equation with little error as 
— Qr, it es h. 
te a 


Finally, if we let C=1Q+C, 
jas =p ewe’ 


and 
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é (fy. Dig 15% Noe? Oily = ee 


icine GE at emmnery cee rags 
or 
hls So) of Waive tee 8 2aryerdNe h,—h- 
, secs aie +2 sf a PO 
of 0 rage: 740,76“ ) 9K h.h_ Q 


With the values chosen this gives a potential difference of —41.2 milli- 
volts, of which —387.5 are across the membrane. 

If we consider a respiratory system such as discussed by A. T. 
Cameron (1930), and H. D. Landahl (19389), 


a g 
glucose + phosphate = hexosediphosphate = lactic acid + phosphate 
(1) (6) (7) (2) (6) 
a 
carbon dioxide + water < oxygen, 
(4) (5) (3) 
we have 
hie ae aCciC, ’ 
Qo nfC.Cs as gC; , 
Ce 1CiGq ; 
C6 =a. QO; = 0, C, + 9C,=0, 
C,= CG; 
We may obtain 
d2— NQ3 — Qi y 
eh N ESF 
Q. = nQ3; A i. Qi. 


If one mol of (1) is oxidized, 38 are glycolyzed and 36 + 1 are used. 
But for oxidation of 1 mol of (1), 60. are used, and 


Q; = 6Q, for oxidation, 


or 
= (3p +1) &, 
and 
_Q:. (88 + 1)Qs 
meer ie roy 
since 
M, 
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Then 


If we identify Q with lactic acid, and assume for a preliminary in- 
vestigation that it is highly ionized and that there is no adsorption of 
either ion, or that the diffusion coefficient is that of the adsorbing par- 
ticle after adsorption, we have then ; 
ae _ AT A. gee kTB A. s a5, 
o OeeteGs 6e Co 


fkTB 


= + = (Ax ad A.) mC; ; 
E 


roid ro SP Ar eAN ch oh. 
P.D.| =PD. (l+—) + ees a Os: 
(0) 0 


9K hh 


2a7rrdNepfmh,—h = = 
RST ea OCs: 

K. J. Lund (1931) gives data relating electrical polarity potentials to 
oxygen consumption at different temperatures. Since his oxygen con- 
sumption is measured as an average over a three hour period, we can- 
not say that the peaks and depressions in the polarity potential curve 
are not due to a variation in O, consumption, since such variations 
were not determined. In fact, we find indicated just such a variation 
with O, consumption if we consider the probable variation of O, con- 
sumption with temperature. For T. C. Barnes (1937) cites the experi- 
ments of T. C. Barnes and T. L. Jahn (1933) in which the previous tem- 
perature variations of water had an effect on Euglenia activity, previ- 
ously cooled water having a stimulating effect. T. C. Barnes ascribes 
this to formation of the trihydrol polymer. F. W. Gray and J. H. 
Cruickshank (1935) found a lag in diamagnetic susceptibility of water 
recently melted, of 20 minutes to the maximum and subsequent sub- 
sidence to a constant value, supposedly a hysterisis lag in polymer 
equilibrium. In E. J. Lund (1931), plate 1, we see a peak in potential 
reached about 20 to 30 minutes after the increase in temperature. 
Correspondingly there is sometimes a significant depression below the 
base line when the temperature is lowered, rather than a monotonic 
approach to equilibrium of the same time interval. This is also ob- 
served in Douglas fir (Lund, 1932a, 1932b). We see then that Lund’s 


= P.D. 


3 %) 
at) - 
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data could be interpreted as a correlation of potential difference and 
metabolic activity in which the lag in activity with temperature varia- 
tion is due to a hysteresis lag in attaining equilibrium of the water 
polymers. This would support our theory, but more information as 
to the variation of O. consumption, P.D., and temperature is needed. 
This effect might also be explained as the result of an accumulation of 
intermediate products at low temperatures which would create an ex- 
cess activity for a brief period after return to a higher temperature. 
However, attempts to store free oxygen and other substances in cells 
have been relatively unsuccessful. 

We have shown that continuously maintained potentials result- 
ing from diffusion are of the correct order of magnitude, and vary in 
the proper way with cell metabolism. Such potentials may play a role 
in nerve excitation, protoplasmic streaming, flow of sap in plants, etc., 
but the forces are too small to be of importance in cell division when 
compared to diffusion forces. These potentials will vanish with death 
of the cell as contrasted to potentials derived from colloid adsorption 
or Donnan membrane permeability (N. Rashevsky, 1935, 1938, 1940). 

The author is indebted to Dr. N. Rashevsky for his frequent and 
invaluable aid, and to Dr. A. S. Householder for checking and correct- 
ing the manuscript. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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NOTE ON THE MATHEMATICAL BIOPHYSICS OF 
TEMPORAL SEQUENCES OF STIMULI 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


Some general considerations are given regarding the effects of 
temporal sequences of stimuli in a neuronice network, which consists 
of a set of parallel chains of excitatory fibers with cross-connections 
made of inhibitory fibers. It is shown that, in general, the excitation 
produced by any individual stimulus of the series is a function of the 
order and duration of the previous stimuli, and that the effect of 
each stimulus thus depends on the whole temporal pattern considered. 


In the mathematical biophysics of the central nervous system 
hitherto only the effects of a number of stimuli applied simultaneous- 
ly, have been considered. The only exception is an important paper 
by H. D. Landahl (1940), in which a special case of two temporally 
separated stimuli are considered. It has been shown elsewhere (Rash- 
evsky, 1938, hereinafter referred to as MB), that in general the prob- 
lem of temporal sequence of stimuli reduces to a problem of a set of 
rather complex integro-differential equations (MB, p. 238). The very 
complexity of those equations and their rather unusual mathematical 
form, makes their usefulness for practical purposes at present rather 


FIGURE 1 
To each fiber of type J corresponds a different stimulus S;. The synapses be- 
tween fibers of type I and fibers of type IJ are marked by s;. In this figure we 
have s, , s, and s,, corresponding to i = 1, 2, 3. A fiber of type II, corresponding 


to the synapse s; is referred in the text as fiber IJ; . 
89 
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doubtful. The purpose of the present paper is to outline a somewhat 
more restricted problem, the mathematical treatment of which does 
not involve any complexities at all, and which at the same time pre- 
sents a definite biophysical interest. 
Consider the usual scheme of a set of parallel excitatory chains, 
with inhibitory cross-connections (Figure 1). The properties of this 
simple structure have been studied from many angles and applied to 
varied cases (MB, chapter xxii ff.). Let us consider at once a more 
general case, by assuming that while the constants of all the excitatory 
fibers are the same, those of the inhibitory fibers are all different. If 
we consider, under these conditions, two stimuli, S, and S; , applied 
to two different fibers of type I, then the inhibitory effect Jia of Sy 
upon the synapse s; corresponding to S; is not the same as the inhibi- 
tory effect Ji, of S; upon the synapse s;, corresponding to S;, even 
when the intensities of S,; and S; are equal. Furthermore let us con- 
sider the case, that all fibers are of the general type, every fiber pro- 
ducing both e and 7. The excitatory fiber produces an excess of «, 
the inhibitory — an excess of 7. Let the former be characterized by 


Av>B 6 307 oe (1) 
while the latter are characterized by 
AP ans. 
A<B; a<b; rire (2) 


Let now a stimulus S; of constant intensity, sufficient to produce 
an excitation in the corresponding fiber JJ (Figure 1), be suddenly 
applied and sustained for some time. This results in an excitation of 
one corresponding center at the efferent end of fiber JZ, and in an 
inhibition of all synapses corresponding to other stimuli S;. Let now 
at some moment ¢ the stimulus S, be suddenly interrupted. With in- 
hibitory fibers characterized by inequalities (2), this will result in 
temporary excesses of « over j , produced by all inhibitory fibers, stim- 
ulated by S;,, at all synapses, corresponding to other stimuli (MB, 
p. 228). The intensity of this “rebound phenomenon” will depend at 
each synapse upon the physical constants of the inhibitory fiber pro- 
ducing it, and upon the intensity of S;,. The time of onset of this re- 
bound phenomenon, at which « becomes greater than 7, will be inde- 
pendent of the intensity of S;,. It will, however, depend on the dura- 
tion of S;. 

Thus the mere cessation of a stimulus results in certain excita- 
tory phenomena, and by an appropriate neural mechanism may be 
perceived itself as a stimulus. 
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For the present let us assume that the rebound phenomenon at all 
the synapses is not strong enough to excite any of the fibers 17. Under 
those conditions the cessation of S; will not excite any of the fibers 
II. But if after its cessation, S; is followed shortly by another stimu- 
lus S; , (applied, of course, to a different fiber of type J) then the ex- 
citation of the fiber IJ, corresponding to S; will be greater, than if S, 
were not preceded by S;,. Inversely, the effect of S; is altered by mak- 
ing it be preceded by S;. Due to Ji: 2 Jix, the two effects are, however, 
not identical quantitatively. The maximum value of excitation pro- 
duced in the two corresponding fibers IJ by the sequence S; S; is dif- 
ferent from the maximum value produced by the sequence S; S;, and 
in general they are both different from the maximum value produced 
by a simultaneous application of S, and S,. Moreover, the effect of 
the combination S; S; will depend both on the duration of S; and S; 
and upon the time interval between them. The rebound phenomenon 
at each synapse starts at zero, reaches a maximum and then drops 
again to zero. Therefore, a maximum increase in the intensity of ex- 
citation of JJ, due to a preceding S; will be reached when a definite 
time t* elapses between the cessation of S, and the initiation of S;. 
This time t* is, however, not identical with the time 4, at which the 
rebound value of « — 7 reaches a maximum (MB, p. 228), because «&; 
and 7; do develop gradually after the application of S;. 

Consider now a sequence of stimuli S; , S;,, S; --- S,, each having 
a corresponding duration ¢;, t&, t --- t,. Denote the time interval 
between S; and S; by ti. Denoting the rebound intensity at s,; due to 
S; by Ry, the quantity Ry: will bea function of S,, % and t. The 
intensity of excitation of JJ, will be determined by all the previous 
S’s, their durations and the intervals between them. For that inten- 
sity is determined by the amount 


Ef = fer Sy ltors (3) 
Dp 


the summation being taken over all previous stimuli. 

A still more complex situation obtains when we consider the 
scheme of Figure 2, in which the inhibitory fibers start at the synapses 
s; (MB, p. 235). In that case even every R,, depends upon all previ- 
ous S’s and all previous R,,’s. We have here an analogy to the phe- 
nomenon, emphasized by many psychologists, namely that the effect 
of a stimulus in a series of others is determined by the “time gestalt” 
of those other stimuli, or, if we may say so, by its “context”. 

In the case represented on Figure 1, the effect of a stimulus is 
always enhanced by preceding stimuli, the amount of enhancement de- 
pending on the past sequences of stimuli. In the case, however, of 
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FIGURE 2 


Figure 2, a preceding stimulus may reduce the response to a given 
stimulus. For the cessation of a preceding stimulus results in a re- 
bound excitation at all synapses, and this rebound excitation excites 
the inhibitory fibers, leading to the synapse under consideration. 
That synapse has now a rebound excitation from the synapse corre- 
sponding to the preceding stimulus and an inhibition due to the re- 
bound excitation at all other synapses. The total “time gestalt” of the 
preceding stimuli will determine which of the two effects prevails, and 
whether the effect of the given stimulus is enhanced or reduced by 
previous stimuli. 
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A NOTE ON THE NATURE OF CORRELATIONS BETWEEN 
DIFFERENT CHARACTERISTICS OF ORGANISMS 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


Different anatomical and physiological characteristics of organ- 
isms affect their interreaction with the inorganic world as well as 
their mutual interreactions. In this way they all may affect indirectly 
the total rate of reproduction of a species. It is shown that the re- 
quirement of a maximum rate of reproduction defines the distribution 
functions of the different characteristics and through those distribu- 
tion functions determines statistical correlations between the char- 
acteristics. 


The relation between different characteristics of an organism is 
in many cases due to the biophysical constants of the organism itself. 
If two different characteristics x and y are determined by the same 
physico-chemical factors, there will be a definite quantitative relation 
between those two characteristics. If, however, the two characteris- 
tics are determined in general by two different groups of factors, 
with only a few factors being common to the two groups, we shall 
find instead of a definite mathematical relation between x and y, only 
a certain correlation, which will be the closer to unity, the more com- 
mon factors the two groups possess. The future development of the 
mathematical biophysics of the organisms as a whole should thus 
eventually lead us to the theoretical prediction of the correlation co- 
efficients between different pairs of characteristics of an organism. 
Such a theory would be of particular value because it will render 
useful a wealth of empirical data accumulated in this field by bio- 
metrists. 

It must be pointed out, however, that there is another possible 
cause for the appearance of correlation between different character- 
istics of an organism. This cause lies in the effects which these given 
characteristics may have upon the interaction of different organisms 
and thus indirectly upon the preservation of the species. Certain pairs 
of values of the two characteristics may be more favorable to the pres- 
ervation of the species. Organisms characterized by such pairs of val- 
ues will have a better chance for survival and such pairs of values will 
therefore occur more frequently than one would expect on the basis 
of a random distribution. The purpose of this note is to discuss briefly 
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the general type of mathematical problems, to which considerations 
of the above mentioned kind may lead. 

Let us first consider the case, that the two characteristics « and 
y vary discontinuously, and may have only discreet values, 7, %, 
Xs, °°' Uns Yr» Yo, Ys? Yn- Let the number of organisms, character- 
ized by the values, xi, y, be equal to Nix. We shall refer to an or- 
ganism, characterized by x;, y, aS organism (ik). The values x; and 
yx, determine amongst other things, the interaction of the particular 
organism with its inorganic surroundings; and through that interac- 
tion they determine the rate of net increase aizN ix, with respect to time 
of the organism of that kind. But that rate of increase also depends 
on the interaction of the organism (7k) with other organisms (7k) as 
well as organisms (m7), where m and n are different from 7 and k. 
We may set that part of the rate of increase as being of the form 


Nix dD OF" Nan. (1) 
We thus have for the total rate of increase of all organisms 
N= 2>2Nix (2) 
+k 
the expression 
dN 
GE TU Ba Nin + DV INin DE OF? Ninn. (3) 


The natural selection acts so as to maximize expression (3), and 
the conditions of the maximum of (3) give us the values of Nix’s. 
These we determine by 


so Tg (SE ee Ne +EINadD be Non) =0 stig 


and this gives a system of equations 


Ain +> > OM" Ninn = 0. (5) 


By a procedure familiar in the theory of integral equations, we may 
pass to the case of continuously varying x’s and y’s. In this case the 
@ix’s and the Nmn’s become functions a(x,y) and N(é,¢) of « and 
y,and and. The coefficients b”” become functions b(a,y,&,¢) of 
four variables. The system of linear equations (5) reduces to the in- 
tegral equation: 


a(x,y) +ffb(z,y,&,0)N(E,¢) dédt=0. (6) 
where a(x,y) and b(x,y,&,f) are known. The solution of (6) 
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gives us the distribution function N(x, y) for both characteristics. 
This function N(«# , y) determines the average values x and y of 
the characteristics. The correlation coefficient between x and “is 
then given by 


(w7— x) (y—y) . (7) 
<i (x — x)?(y — y)? 
and is completely determined by N(z, y). 

The solution of the integral equation (6) represents the mathe- 
matical problem involved in the development of the suggested theory. 
The biophysical problem lies in the determination of a(x,y) and of 
b(x,y,&,¢), and represents a different problem for every different 
meaning of x and y. For instance, x may stand for the weight of the 
animal, and y for its speed of locomotion. In determining b(x,y,&,¢) 
we will have to consider that a large x is beneficial for the preserva- 
tion of the species by offering a protection against smaller enemies. 
A large y is also beneficial, by enabling avoidance of enemies on one 
hand and catching prey on the other. But the effect of € and ¢ is op- 
posite. In general in this case b(x,y,&,¢) will increase with x and 
y and decrease with € and ¢. 

The salient feature of the situation is that the characteristics of 
an organism depend on the structure of the whole organic world. A 
particular distribution function N(z,y) is determined by its con- 
tribution to the maximalization of the rate of reproduction of the 
organic world as a whole, and not of a particular species. We thus 
may find that certain characteristics of an organism exist, which are 
of no apparent use to that organism (C.f. appendix and caecum). 

We may set approximately b(x,y,&,¢) and a(x,y) as linear 
functions of x, y, & and ¢, with indeterminate coefficients, and then 
determine these coefficients from the observed correlations. 


T= 
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Beginning with Rashevsky’s equation for the development of 
the excitatory state in a nerve fiber, an equation for the change in 
skin resistance upon the presentation of an instantaneous stimulus is 
derived. The mechanism assumed is in conformity with the existing 
evidence of neuro-physiology. Certain deductions from the equations 
are made and experimental problems suggested for testing the theory. 


One of the most perplexing problems faced by the physiological 
psychologist is the problem of how to measure the galvanic skin re- 
sponse. Upon presentation of a stimulus, a curve indicating the change 
in the apparent skin resistance is secured. The basis of the problem 
lies primarily in the fact that there is no clear idea of the nature of 
the biophysical mechanism mediating the response and hence no clear 
idea of the significance of the various parameters of the curve. 

There is, however, sufficient evidence (Richter, 1929a; O’Leary, 
1932; Darrow, 1934; Dale, 1986; Rashevsky, 1938) from a variety of 
sources which may be integrated to provide a reasonable hypothesis 
for the biophysical mechanism. Deriving an equation for the change 
in skin resistance from a quantitative expression of the hypothesis 
will, when fitted to the data, yield parameters that have physiological 
and psychological significance. 

That the galvanic skin response undoubtedly has a very complex 
mechanism with a variety of factors interacting is agreed upon by all 
workers in the field. But it is not until quantitative expression is 
given to hypotheses as to the nature of the interaction that the valid- 
ity of the hypotheses may be checked against experimental data. 

In a first attempt to develop a quantitative theory, it would be 
futile to try to include all of the great variety of factors that con- 
ceivably might enter into the determination of the finally observed 
response. In a first approximation we shall select only those factors 
which we consider to be of fundamental importance and impose the 
simplest experimental conditions upon the equations. If the experi- 
mental data are in general agreement with the predictions of the 


97 


98 MATHEMATICAL BIOPHYSICS 


equations we may then feel that we are on the right track. The next 
step is to complicate the experimental conditions, making no change in 
the postulates in regard to the biophysical mechanism, until a point 
is reached at which the experimental data fails to conform to the 
theory. This would then indicate that one or more of the factors pre- 
viously neglected in the development of the equations has, with the 
more complex experimental conditions, become important; and a sec- 
ond approximation or extension of the theory is in order. 

On the basis of the sources previously referred to, the biophysical 
mechanism hypothesized as underlying the galvanic skin response is, 
in its simplest form, as follows. Upon presentation of a stimulus 
above some threshold an excitation of the autonomic nerves mediat- 
ing the galvanic skin response occurs. The intensity of the stimulus 
applied to the autonomic nerves we shall designate by the letter EF. 
An excitatory state (e) is developed in the autonomic nerves at a rate 
proportional to FE’, and is dissipated at a rate proportional to itself. 
A chemical mediator (p,) is produced at the peripheral end of the 
autonomic fibers as a result of the excitatory state, and diffuses to the 
sweat gland tissue, in the neighborhood of which it is designated as 
pa. We shall further assume that pz bears an inverse relation to the 
apparent skin resistance. 

This hypothesis may be developed quantitatively in the following 
manner. 

The development of ¢ is given by N. Rashevsky (1938, chap. xxii) 


de : 
Ts =A'H—ae. (1) 
Introducing a change of variable 
at=-r 
A A’ 
a 
adt = d i 
and substituting in (1) 
de 
re =AEH-e. (2) 
The quantity « varies according to 
e=AEH(1—e-), (3) 
In the absence of any stimulus EZ , equation (2) takes the form 
de 


ieee ie (4) 
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and the quantity « varies according to 
e=6,e7 (5) 


where e, is the value which « has at the moment the stimulus ceases. 
The simplest experimental condition which we may impose is 
the case of an instantaneous stimulus, where the duration of the stim- 
ulus will be 7,. Although 7, is negligible, H 7, must be finite for a 
change in skin resistance to occur. 
The initial condition ¢, is the value of « in equation (3) attr = 7, 


o — AP (1 — 6-73), 
Expanding (1 — e”) and taking only the linear terms in 7, leads to 
é,—AE7,=¢q, 
and this substituted in (5) gives 
e=qe7 (6) 


which is the value of ¢ for 7 > 7,. 

Since 7, is small, the chemical mediator, p,, will remain un- 
changed for 0 <z<7,. But for 7 > 7, it develops at a rate propor- 
tional to « and is dissipated at a rate proportional to itself. Hence, 
we have 

ce — fb bps (7) 
es 
where é and wp are constants of proportionality. 

Substituting (6) in (7) and dropping all terms in 7, of higher 

order than the first, we have for the value of p, 


Pp = = i (es (8) 
b 


This is the value of the concentration of the chemical mediator 
of excitation at the peripheral end of the nerve. 

- The chemical diffuses to the region of the sweat glands where it 
modifies their permeability directly in proportion to its concentration. 
Designating its concentration (or the permeability of the sweat 
glands) as pz, and letting p, represent the inherent permeability of 
the tissue in the absence of a chemical mediator, we have 


d (Da — Po) he 


2 O Dp — (Da — Po) (9) 


where a and » are constants of proportionality. Inasmuch as pp ‘and 
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(pa — Po) represent the same substance or state, whatever it may be, 
at two different places, we have, as a first approximation, assumed 
their constants of proportionality (w) in equation (7) and (9) to be 
equal. 

Substituting (8) (9), integrating and defining 


P= patr=71 


where p designates residual permeability remaining from previous 
stimulation, we obtain, after dropping all terms in 7, of higher order 
than the first, 


pi= po + (p— poet + 4 [er — [1- (1—w)r]e*|. 10) 
(ESky 3 
If we assume that conductance (c) bears a direct linear relation 
to pa 
pe 
k 


we obtain, by substituting in (10) and rearranging, 


eae! ete ea eee em. (11) 


The equation describing the variation of resistance (7) with 
time may be easily secured by substituting ¢c = 1/r. 
It is readily observed that at 7 = 0 


Pp 
Co Sarit (12) 
and as +r > 
ox Be (13) 
k ? 
while the initial slope 
o) =-£ @ =p) =— mle — en) 38. (14) 


From a consideration of equations (12)-(14) certain interesting 
experiments suggest themselves. 

1. Suppose a succession of instantaneous stimuli are presented 
the subject. For stimulus (i) the initial slope (S;) and the initial 
conductance (c¢.;) may be obtained. Substituting these values in equa- 
tion (14) we have : ae 
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SC, = Coo) =S8; , 
and similarly for stimulus (7) 

= BCG Con = S 
Taking the ratios and solving for c., we obtain 


Coj9i 7 CoiS; 


35. (15) 


es eee, 
As a consequence of the assumptions so far introduced, ¢.,, ob- 
tained by equation (15) from the responses to any two stimuli (i, 7), 
should be constant for a single subject. However, there is reason to 
believe (Richter, 1926b, 1929b), especially in the waking subject, that 
there is some degree of continuous excitation of the autonomic fibers 
mediating the galvanic skin response. Hence the excitation induced 
by stimulation of the subject in the laboratory will be superimposed 
on the more or less constant excitation existing at the moment. 
Thus as 7 > o 


Po , D 
Cae ke 
where ’ represents the state of permeability of the tissue induced by 
the more or-less constant level of excitation that prevails in the sub- 
ject. 

An alternative method of introducing the concept of p’ is to in- 
troduce a positive constant E’ in equation (2) and carry through the 
integrations in the same manner as before. Equation (15) when ap- 
plied experimentally will yield a constant only if p’ is constant. If p’ 
varies in the course of an experiment, the values obtained in equation 
(15) will vary directly with it. 

It may even be suggested that the lability of the parameter p’ 
may be related to certain temperamental characteristics of the subject. 

2. In order to obtain some qualitative estimate of the general 
nature of the time relations and the parameters of the function, we 
may differentiate equation (11) with respect to 7, and set it equal to 
zero. Doing this and rearranging terms, we obtain 


(Dp — Do) (u — 1)? 


aégq Ce 


e pt — 4(1 — a) rc + 1 


If we transform the time-units by 


e£=(1-4p)r, 
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we have 
D — Do) ( oe 
ae 2 Pate ie (17) 
aéq 
Then if we set 
ye", (18) 
Dp — ay ¢ cee Fe 


acd 

the intersection of the exponential curve and the straight line gives 
the value of x at which we have a maximum or minimum in the gal- 
vanic skin response. The following relations may be seen to hold in 
equations (18) and (19). For p > p., other things being equal, as 
(p — po) increases, the latent time increases and the time to maximum 
conductance decreases. An approximate expression for the latent time 
under these conditions may be obtained by taking the first two terms 
of the expansion of the left hand side of equation (17), solving for x 
and substituting (1 — w)r— «x. Doing this, we have 


Ey ps) 
—— Wer 
where 7’ is the latent time. 
For p < p., the latent time cannot be secured from the equations, 
but as p decreases the time to maximum conductance increases. 
Also, it may easily be seen from equations (18) and (19) that, 
other things remaining constant, for p > p., aS q = AE 7, increases 


the latent time decreases, and the time to maximum conductance in- 
creases. 


For p < p,, aS ¢g = AE7, increases the time to maximum conduc- 
tance decreases. Nothing may be said with regard to the latent time 
for p< p,. 


For p = p,, the time to maximum conductance is independent of 
the stimulus intensity. 


3. For p> p,, if a subject be given successive equal stimuli (2), 
such that each stimulus occurs a few seconds after the point of maxi- 
mum conductance is reached, but before recovery is complete, the 
initial slope (coi) of each successive response will be numerically 
greater, the magnitude of the response, measured in conductance, less 


than the preceding, the latent time will increase, and the time to 
maximum decrease. 
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From equation (14) it is obvious that as (p — p,) increases, the 
initial slope increases numerically. 

That the latency will increase and the time to maximum decrease 
has already been shown. This indicates that the duration of the ex- 
cursion decreases, and it is reasonable that its magnitude also de- 
creases. The limiting value for the duration of the excursion is zero 


as (p — p.) increases, hence the limiting value of the magnitude of 
the excursion is zero. 

For a given q, if p is very large, neither minimum nor maximum 
can occur. This means that the conductance can only decrease, despite 
the applied stimulus. Increasing the stimulus, however, will lead again 
to a maximum. 

Inasmuch as gq represents essentially the intensity of the stimu- 
lus applied to the autonomic fibers, it may or may not be a direct cor- 
relate of the physical intensity of the stimulus. With frequent suc- 
cessive presentations of the same physical stimulus, there may be some 
sort of central adaptation decreasing the effective value of g (Coombs, 
1938). 

The author wishes to acknowledge his appreciation of the aid of 
Dr. N. Rashevsky, who made this study possible, and of the many sug- 
gestions and criticisms of Dr. Alston S. Householder. 

This work was aided in part by a grant from the Dr. Wallace C. 
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It is found that for a simple circuit of neurons, if this contains 
an odd number of inhibitory fibers, or none at all, or if the product of 
the activity parameters is less than unity, then the stimulus pattern 
always determines uniquely the steady-state activity. For circuits 
not of one of these types, it is possible to classify exclusively and 
exhaustively all possible activity patterns into three types, here 
called “odd”, “even”, and “mixed”. For any pattern of odd type and 
any pattern of even type there always exists a stimulus pattern con- 
sistent with both, but in no other way can such an association of ac- 
tivity patterns be made. 


In continuation of an investigation previously introduced (House- 
holder, 1941; this paper will be designated “I’’, hereafter), we con- 
sider here n simple neurons (each having a single axon and a single 
dendrite) which we shall designate (1), (2), -:-, (m), arranged in a 
simple circuit so that fiber (7) synapses with fiber (¢ + 1). We shall, 
when convenient, use the notation (7) and (n + 7) interchangeably to 
denote the same fiber. We suppose that if a constant stimulus S; is 
applied to neuron (2), and if S; exceeds the threshold h; of this neu- 
ron, then, in the absence of other factors, this neuron tends asympto- 
tically to a steady-state activity in which it develops « — 7 in the 
amount a;(S; — hi). (For terms and symbols not here defined see 
Rashevsky, 1940, and I.) 

It should be pointed out that whereas we speak here of single 
neurons reaching a steady-state of activity, this is not the only pos- 
sible physical realization of the equations and relations assumed. Thus 
it may be that a small bundle of several parallel fibers is involved, in- 
stead of a single one, so associated that while adaptation leads to a 
periodic inactivity of one or another fiber of this bundle, the activity 
is taken over by other fibers of the bundle, the net effect being one of 
steady activity. Again, an inhibitory fiber may be a fiber or bundle of 
fibers so related to the next fiber or bundle that its activity serves 
always to depress these. In this connection cf. N. Rashevsky (1940, 
chap. vili). 

, The constant a;, which may be positive or negative, but is not 
null (when any a; = 0, there is no problem), is called the activity para- 
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meter of the neuron. Hence, if at each synapse a constant stimulus 
S; is applied from a source which is independent of the activity of 
neurons of the circuit the total stimulus 7; acting upon (7) is the sum 
of S; and the contribution of the fiber (¢ — 1). By setting 


Yi=ni—h, o=—S;—h,, 
as in I, the activity of the circuit is described by the equations 


Yi =o; T+ Gin Yia» (1) 
where 
aj =—0 when %=0, 
(2) 
a;= a; when. 4; > 0. 
The vectors 


= (Gy ,°++, On), Y = (Yi5°++5 Yn), o = (01,°**5 on) (3) 


we shall call, respectively, the activity pattern (AP), the excitation 
pattern (EP), and the stimulus pattern (SP) of the circuit. This 
terminology deviates slightly from that employed in I. Evidently y 
determines uniquely a and o; also o and a together determine uniquely 
y. But not every pair of vectors « and a determine a y by (1) which 
is consistent with (2). When they do, we shall say that o and a are 
consistent. In fact, one result in I was to the effect that for circuits 
containing only excitatory neurons (all a; > 0), there is one and only 
one a consistent with any given o. 

It is evident from (2) that a is known when its null components 
are known, since its other components are equal to the corresponding 
a;. If o, = 0, we shall say that a, or the AP, has a zero at k, or that 
k is a zero of a. Another result of I was to the effect that if two vec- 
tors a and f are AP’s having a common zero, and if both are consis- 
tent with some oa, then a = £. We shall consider here only AP’s a with 
at least one zero, and the problem is, for a given circuit, if a has at 
least one zero, how many AP’s £ are there, each of which has at least 
one zero, and each of which is consistent with some o with which a 
is itself consistent? Or, given o, with how many AP’s, each contain- 
ing at least one zero, is it consistent? 

Thus along with (1) and (2), we consider 


204 + Bin 2a (4) 
where 


8; =0 when 2,50, 


Bj; =a; when 2;>0. (5) 
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We suppose a # 6. Hence no k is a zero of both ao and 6. As in I, we 
define 


a = Tis 
o =o + ay, 0), (6) 
7 7 4-1 
Hence 
0) =o!) + Oia Gig + Op 0). (7) 
From the matrix (o™) pick out the set 
» x An 
of ae ae see ) (8) 
by means of the conditions 
Ais =1-+ hi when Qj; #0 ; 
(9) 


crea | when oa,=—0. 


Then, as we have seen in IJ, if a is, in fact, consistent with o, it must 
be true that 
7.00 (10) 


satisfy (1) consistently with (2), and hence 


ae) = 0 when ,a;—0, 


ao?) >0 when a,—4a;. GD 
In like manner we determine the set 
61 P2), GaP) , ++, y(t (12) 
by means of the conditions 
Mins = 1+y; when £40, 
Hin = 1 when §,=—0; ee 
equations (1) have the solutions 
2; = oH (14) 
consistent with (5), and 
oi4) =0 when 6;=—0, 
(15) 


oi") >0 when 6,40. 


The operations just carried out may be conveniently represented 
as follows. Consider a cylindrical surface of circumference n and al- 
titude n, and the lattice of points with positive integral coordinates. 
The lattice point whose abscissa is i and whose ordinate is 1 may be 
associated with the element o; of the matrix of o’s. With any AP, 
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a, we may associate a curve on this cylinder, the curve consisting en- 
tirely of helical arcs joining the points (i, 4;) as the 4; are defined by 
(9). Thus all ascending arcs of the curve become lines of unit slope 
when the cylindrical surface is slit along a generator and flattened 
out onto a plane. And every arc is an ascending arc except those lead- 
ing from points of abscissa 7 for which a; = 0. Hence if a; = 0 we 
may speak of the point (i, 4:) as being a (relative) peak, and it is an 
actual peak unless also a;-. = 0. Thus every relative peak of this 
curve marks a zero of the AP, a, and conversely. 

Associated with any other AP, 8, will be a similar curve, and 
these two curves may be called the a-curve and the f-curve respective- 
ly. Then if o and § are both consistent with some oa, no peak of the 
a-curve can have the same abscissa as a peak of the f-curve. Hence 
while the two curves must cross, they can have no common arc and no 
common lattice point. This representation is not essential to the en- 
suing discussion, but it may serve to clarify some of the statements. 

We prove first 


LEMMA 1. If the AP’s a and §, each containing at least one zero, 
are both consistent with the SP o, and if a (or B) has zeroes at h and 
k , then 

Qn Unser °°* Ana = Og A Ans *°* Ua >0 ’ ; (16) 


while if h and k are zeroes, one of a, the other of 8, then 
Qn Aner? Ona <0, A, Aer °° Ma <0, (17) 
any subscript exceeding n being reduced by n. As a special case of 
(16) it follows that 
Cilgk a Ogi 0' 2 | (18) 


First, suppose 7 and 7 + y are zeroes of a, but that the sequence 
a+ 1,7+ 2,---,7%+»— 1 contains no zero of either a or £. Then it 
follows from (9) and (13) that 


disp =p» [isp = Mi + p (p=1,---, 7); (19) 
from (11) it follows that 
ws , a7 0 (p=1,---,»—1) ’ (20) 


and from (15) that 
mee >0 (p=1,---,¥). (21) 
Finally we have from (7) that 


(pitP) —= (p) F : see Oe = 
Finp Tip Qi+p-1 Qisp-2 a; ok (p = i Ss y) . (22) 
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With h =i, k =i + », the first of (16) follows in this. special case 
from (20), (21) and (22) when p = ». 

Next suppose that i + » + »' is the next zero of a, and-that the 
Setity+1,it+v4+2,---,i+»+»' — 1 contains no zero of either 
aor 6. Then 

Q; Vin? Misa > 0 » Visy Wisya s+? eae 1 > 0 ’ 
and therefore the first of (16) follows withh =i,k=i+»+». By 
a simple induction the first of (16) can be proved in all cases where 
the seth +1,h+2,---,k—1 contains no zero of £. 

Before proceeding with (16) we consider (17). Let i be a zero 
of a, 2+ » of 6, and suppose that the sequence i + 1,7 + 2, 

7 + » — 1 contains no zero of either a or 6. Then (19) and (22) still 
hold, but in place of (20) and (21) we have 


ea (p=1,---,”), (23) 
reat =O, Pea >0 (p=1,---,»—1). (24) 


Then with h =i, k =i + », the first of (17) follows in this special 
case by taking p = » in (22), (23) and (24). But now the induction 
can be completed, for if we pick out from the sequenceh,h +1, 
k—1,k, all those terms h, =h, h.,--+, hy =h, each of which is a 
zero of either a or f, then in the hypothesis for (16) there will be an 
even number of pairs of consecutive terms of this subsequence, h, and 
hy. , of which one is a zero of a and the other of 6; while under the 
hypothesis for (17) there will be an odd number of such. 

COROLLARY 1. If the circuit contains an odd number of inhibi- 
tory fibers (the number of negative a; is odd) or if it has none, then 
for any SP o there cannot be more than one AP a which contains at 
least one zero and is consistent with o. Moreover y and a can be deter- 
mined from the matrix of the o,. 

We, therefore, consider only circuits which contain an even num- 
ber of inhibitory fibers in the sequel. Let there be 2m inhibitory fibers 
and u = n — 2m excitatory fibers. Suppose the Beate fibers are 


(41), (%2), cea: (tom) 
and the excitatory fibers 


| (31), (ia)o-e > Cu), 
and suppose that 
2 nating LSi<i<-+ <twS 0, ' 
re19 Aes Galiano 3p kG Be ath oa 
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Let a be any AP with at least one zero and let k be one of its zeroes. 
Then there is a smallest integer «x (which may be zero) for which 
k+.«(ork+x«—7) is equal to somei,(r=1,---, 2m). We shall say 
that the zero at k, or that the fiber (k) is of even type when r is an 
even number, and that it is of odd type when 7 is an odd number. Evi- 
dently a change in the enumeration of the neurons will either inter- 
change all types or else leave them all fixed. If a has zeroes of even 
type only then a will be said to be of even type; if a has zeroes of odd 
type only then « will be said to be of odd type; if a has zeroes of both 
types, a will be said to be of mixed type. Then it is only a restatement 
of lemma 1 to say 


LEMMA 2. If the AP’s a and B, each containing at least one zero, 
are both consistent with some SP «a, then either a ts of odd type and B 
of even type, or vice versa. In particular, no SP o can be consistent 
with more than two AP’s and if it is consistent with any AP of mixed 
type, then it is consistent with this alone. 

We have yet to see whether, given an AP a of odd type and an 
AP £ of even type, it is always possible to find a « consistent with both. 
Before doing this we must examine further the conditions which we 
imposed upon the o’s when this is to hold. These relations can be sum- 
marized thus: 


LEMMA 3. If a and Bf are two AP’s, each containing at least one 
zero and each consistent with the SP o, and 


a) tf hand kare zeroes of a, then 
Ze = An U2 +++ Uy Zn} (25) 


b) «if hand k are zeroes of 8, then 


Yur = Ua ++ Un Yr; (26) 
c) tfhisazero of a and k is a zero of B, then 
Yr. = Apa -++ Uy Zn5 (27) 
d) ifhisa zero of B and : a zero of a, then 
Ze S — Anas Yn. (28) 


The y; and 2; are understood to be the solutions of (1) and (4). 
Bearing in mind equations (10) and (14), we proceed as in the 
proof of lemma 1. If 7 and i + » are zeroes of a, but the sequence, 
1+1,---,t +» — 1 contains no zero of a or of 6, then (19) —(22) 
hold, and these, for p = », give (25) in this special case. If i is a zero 
of a,i+ » of 8, and the sequencei + 1,---,i + »— 1 contains no zero 
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of either a or of 6, then (19) and (22)—(24) hold and we have (27) 
in this special case. By interchanging the roles of a and B we obtain 
(26) and (28) in the corresponding special cases. In the general case 
we pass from zero to zero applying the special formulas, and so ob- 
tain the desired result by induction. 

If we take h = k in either (25) or (26), then since z, > 0 or 
Yn > 0 as the case may be, we have 

COROLLARY 2. If there exist two AP’s, a and 6, each having at 
least one zero, and an SP o consistent with both, then 


aoe. 8 0 12 (29) 


These results can now be summarized and completed. 

THEOREM. Let a circuit be given which has in it at least two in- 
hibitory neurons, and let the activity parameters of the neurons sat- 
isfy (29). Let a be any AP of odd type and 6B any AP of even type. 
Then there always exists an AP o consistent with both of them. But 
for no other type of circuit, and for no other kind of pairing of AP’s 
can such ao be found. 

Having given the zeroes of each AP, we can determine the super- 
scripts 2; and uw; of (10) and (14). At all the zeroes h of a and 6 
there are relations of the form (25)—(28) restricting the values of 
the y;, and z,. These are consistent provided only (29) holds, although 
if the equality sign holds in (29) it must also be taken in (25)— (28) 
and then the choice of a single y; or z, prescribes the values of all the 


others. 
Between pairs of zeroes, 7 andi + », we have relations of the form 


Zisp = Yisp + Vispa es i 2% (30) 
if 7 is a zero of a, and of the form 
Yisp = Zisp + Vispa eG Yi (31) 


if 7 is a zero of f, as follows from (22). By setting p = », it follows 
that having chosen the y, and z, occurring in (25)—(28), the y, and 
z, at the zeroes are also fixed. Now inequalities (20), (21), (23), 
(24) and their analogues with a and f interchanged, and with p # 
when expressed in terms of the y’s and z’s, impose lower limits on 
these, but not upper limits. Hence these are always consistent and for 
each p, either Yi+p OF Zi+p can be chosen, and the other one is fixed. The 
o; are uniquely determined once the y; (or zi) have been fixed. 
We conclude by noting that if, among the pu excitatory fibers (j1), 
--, (jn), there are uw, of odd type and 2 = u — pu, of even type, then 


there are 
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PAwigy o1 | 
distinct AP’s a of odd type and 
D+ 2 ios 1 


distinct AP’s 6 of even type. Hence if (29) holds there are 
(2m — 1) (202 — 1) 


ways in which a pair of AP’s a and # can be found both of which are 
consistent with some SP oc. 

As a simple illustration of the foregoing, consider a circuit of 
four fibers. If the number of inhibitory fibers in this circuit is 0, 1, or 
3, or if a,d.a,0, < 1, then any o is consistent with at most one AP hav- 
ing at least one zero. Suppose, however, that 


Oe> 0 - G2 04--0e:<.0 5) Fe 04 OC 00e ek, 


Then (4) alone is of even type, (1), (2) and (8) are of odd type. 
Hence , = 2, # = 0, y= 1, and there are 2? — 1=7 AP’s of odd 
type, 21 — 1 = 1 AP’s of even type. The latter is of the form (a, @, 
a; , 0); the others are of the form (0,, o, a3, @4), Where at least one 
a is zero, and any one of these is consistent with some o which is con- 
sistent with the single AP of even type, but not with any o which is 
consistent with any other AP. Moreover no AP of the type (a,, ®, 
a3, 0), where at least one of the a’s is zero, is consistent with any o 
which is consistent with any other AP. 

In the circuit for which a, > 0, a.<0,a,; > 0, a, < 0, there are 
three AP’s of odd type of the form (a,, a, a3, @4), and three of even 
type of the form (@,, a, a3, a4), and if any AP of odd type is asso- 
ciated with any AP of even type, a o can be found which is consistent 
with both, provided a, a a; @ 21. 
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